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A systematic analysis of low temperature magnetic phase diagrams of Ce compounds is performed 

in order to recognize the thermodynamic conditions to be fulfilled by those systems to reach a 
quantum critical regime and, alternatively, to identify other kinds of low temperature behaviors. 
Based on specific heat (Cm) and entropy (Sm) results, three different types of phase diagrams arc 
recognized: i) with the entropy involved into the ordered phase (Smo) decreasing proportionally to 
the ordering temperature (Tmo), ii) those showing a transference of degrees of freedom from the 
ordered phase to a non-magnetic component, with their Cm{TMo) jump (ACm,) vanishing at finite 
temperature, and iii) those ending in a critical point at finite temperature because their AC^ do 
not decrease with Tmo producing an entropy accumulation at low temperature. 

Only those systems belonging to the first case, i.e. with Smo — >■ as Tmo 0, can be regarded 
as candidates for quantum critical behavior. Their magnetic phase boundaries deviate from the 
classical negative curvature below T w 2.5 K, denouncing frequent misleading extrapolations down 
to T = 0. Different characteristic concentrations are recognized and analyzed for Ce-ligand alloyed 
systems. Particularly, a prc-critical region is identified, where the nature of the magnetic transition 
undergoes significant modifications, with its dCm/dT discontirmity strongly affected by magnetic 
field and showing an increasing remnant entropy at T — >■ 0. Physical constraints arising from the 
third law at T — >■ are discussed and recognized from experimental results. 
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I. INTRODUCTION 

Magnetic phase diagrams are drown resuming the relevant characteristics of magnetic systems once their basic 
thermodynamic properties are recognized. Their comparison allows to distinguish between general physical phenomena 
and the particular behavior of a single compound. The respective magnetic phase boundaries can be traced by applying 
standard external control parameters like chemical composition (x), pressure (p) or magnetic field (B), able to drive 
the system into regions of relevant interest. Among them, those related with magnetic instabilities, frustrations, 
critical points, exotic phases became increasingly attractive because they involve novel phenomena, stressing the 
understanding of basic physical concepts. 

At low temperature, where thermal and quantum fluctuations compete in energy, the neighborhood of quantum 
critical points (QCP) allows to confront the well established thermodynamical laws with new experimental evidences 
of quantum effects. In fact, a QCP is currently defined as the T = limit for a second order transition driven 
by one of the mentioned non-thermal control parameters [1 . Despite of its unattainable nature at T = 0, a QCP 
carries a sort of halo produced by quantum fluctuations whose physical effects are observed at finite temperature. 
The phenomenology arising from those low lying energy excitations is usually identified as 'non-Fermi-liquid' (NFL) 
behavior |2 , in contrast to the canonical Fermi-liquid (FL) observed in non-magnetic systems. One of the most 
relevant features of NFL systems is the increasing density of low energy excitations, manifested as a divergence of 
thermodynamic parameters like specific heat and thermal expansion divided temperature, and magnetic susceptibility 
when T — >■ 0. Accordingly, also electrical resistivity deviates from the dependence of a FL showing a typical 
linear thermal dependence 0. Within this region, the classical magnetic phase transitions (dominated by thermal 
fluctuations) transform into quantum phase transitions (QPT) gradually dominated by quantum fluctuations [3]. 

In this work, a comparative analysis of the thermodynamic behavior of Ce-lattice exemplary compounds is carried 
in order to determine the conditions upon which quantum critical or alternative behaviors can be expected. After 
to go over the magnetic phase transitions induced by pressure in metallic Ce, an early encompassing phase diagram 
performed on Ce-binary compounds is reviewed. This phase diagram allowed to recognize some relevant concentrations 
for Ce-ligand alloyed systems. The second section is devoted to identify different types of entropy evolutions as the 
ordering temperature decreases because in real systems not all phase boundaries can be extrapolated to zero without 
basic thermodynamical principles violation. The distinct properties of systems fulfilling the conditions to access to 
a quantum critical region are presented and discussed in section III, and the thermodynamic implications of the 
third law approaching the T — !■ limit are analyzed in section IV. Hereafter, the parameter Tmo refers to the order 
temperature independently of its antiferro (AF) or ferromagnetic (FM) character. However, if any discussed property 
is only present in AF systems, the phase transition will be identified by the usual Neel temperature Tpf. 



A. Ce metal phase diagram 

Despite of the significant progress done during the last decades in the study of low temperature phase diagrams 
of intermetallic compounds containing Ce, Yb and U lattices |2], some basic phenomena discovered time ago are still 
under discussion. To our knowledge, the first evidence of anomalous behavior of the 4/-electrons was observed nearly 
eighty years ago when the magnetic behavior of CeN was investigated in 1936 1]. Magnetic measurements showed a 
weak temperature dependence with a fractional value of the Ce-4/ magnetic moment respect to the expected from 
the Hund's rule J — 5/2 angular moment. This observation suggested the concept of "intermediate valence" (IV) 
for the first time in '/' elements. Later, on the early sixties, the phase diagram of Ce metal became the subject of a 
systematic study [S]. Driven by pressure, it displays a diversity of phases where structural and magnetic changes are 
strongly related. The "f -H- a structural and magnetic phase transition became the fingerprint of the local-itinerant 
dilemma of Ce-4/^ electrons not yet completely elucidated. 

That first order transition is related to the collapse of the Ce atomic volume (about 15%) and shows an end critical 
point (CP) at w 600 K under a pressure of ss 2 GPa. Within the a structure, a superconductive phase develops up 
to Tsc ~ 50 mK, which jumps up to Tgc = 1.9 K in a second structural collapse between the a -O- a' phases at 4GPa 
[7] . This superconductive 4/ mediated phase compares in temperature and IV character with the recently highlighted 
second superconductive dome of CeCu2Si2 tuned by pressure and claimed to be related to a second QCP in that 
compound [8]. Since it occurs at the edge of the heavy fermion HF-IV crossover, where the N=2 degeneracy of the 
HF-ground state (GS) transforms into a N=6 GS, this critical region can be regarded as a reminiscence of the a ^ a' 
transition. Including the CeTIns compounds [S' and ZrCe alloys |10j . « 2 K seems to be an upper limit for Tgc already 
detected four decades ago. 
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FIG. 1. Schematic magnetic phase diagrams after [13]. (a) Tmo ordering temperature, xq concentration of the Tmo{x) — >■ 
extrapolation, xk onset of Tk{x) oc 9p increase, LM phase boundary of Ce compounds with local magnetism (xq < xk) 
and ItM with itinerant magnetic character (x > xk)- (b) Encompassing phase diagram, Acf crystal field splitting, Tgs,cf 
respective ground and excited levels hybridization strengths, N different GS degeneracies. 



B. Low temperature phase diagrams of Ce- lattice compounds 

A turning-point in the description of Ce, Yb and U magnetic phase diagrams occurred when the competition 
between on-site Kondo screening and inter-site RKKY interactions was taken into account by theory |llj . While the 
former effect weakens the local Ce — 4/ moments, the latter provides the long range interaction favoring the magnetic 
order. The characteristic energies of these mechanisms {UbTk and ksTR) can be compared through their respective 
dependencies on the local-conduction band coupling parameter (g — SpJex)' oc exp{l/g) and Tji oc g^, where Sp 
is the density of states and Jex{< 0) the usual exchange integral. The nature of the GS is therefore established by the 
value of g, being magnetic (Tr > Tk) for small g and non-magnetic (Tr- > Tr) for large g values. At the intermediate 
region, where Tr- and Tf? are comparable, both mechanisms compete, weakening the local 4/ effective moment with 
the consequent reduction of the ordering temperature Tmo (g) ■ 

Although this model contains the basic mechanisms which govern the magnetic phase boundaries variation, it is 
evident that a 'one-parameter' description cannot cover the diversity of behaviors observed applying different control 
parameters. As example one can mention the different behavior induced by the three usual control parameters: x, p 
and H, on the specific heat of a single sample hke CeCu2(Sio.9Geo.i)2 presente in Fig. 15 of Ref.|T2]. Furthermore, one 
may observe different 'trajectories' of the magnetic phase boundaries between the ferromagnetic CePd (Tc = 6.5 K) 
and two IV isotypic compounds CeRh and CeNi, in both cases driven by alloying the Ce-ligand atom. In the former 
case: Ce(Pdi_3;Rh2,), there is a modification of the chemical potential whereas in the later: Ce(Pdi_a;Ni2.), a structural 
pressure effect [IS . Notice the different effect produced by both control parameters which are usually confused as 
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FIG. 2. Schematic representation of three possible dependencies of Smo{Tmo)'- 1) (sohd hne) fulfills the condition to reach a 
QCP because Smo — >■ as Tmo — ^ 0; 2) (dash-dot line) a phase boundary vanishing at finite temperature and 3) (dashed line) 
a phase boundary ending at a finite temperature critical point (CP), see the text. 

being equivalent [l4] . 

C. Characteristic concentrations in Ce-ligand alloyed systems 

Another limitation for the description of the magnetic phase diagrams using a single parameter dependence arose 
from the comparative analysis performed on seventeen alloyed Ce-binaries compounds [13] • For such analysis, two 
characteristic concentrations were defined: xq, where the magnetic phase boundary Tmo{x) extrapolates to T = 0, 
and xk where the paramagnetic temperature Op starts to rise significantly. It is known that Op oc Tr- [15j once 
|0p(a;)| becomes much larger than Tmo- Two distinct types of phase diagrams were identified, one with xq > xk 
and the other with xq < xk , which are schematically represented in Fig. [T^. Both scenarios correspond to different 
hybridization strengths (F cx Tr-): i) with weak hybridization, i.e. of local moment (LM) character, and ii) with 
moderate hybridization, i.e. of itinerant (ItM) character [16]. This scheme applies to a doublet (N=2 degenerated) 
ground state of Ce-4/ ions within the region where Tmo{x) decreases. For x > xx, T on Op increases more rapidly 
with the consequent broadening of the magnetic levels. Eventually, an overlap between the ground Tgs and crystal 
field (CF) excited levels may occur once respective Tgs a-nd Tcf hybridizations strengths become comparable to the 
CF splitting (Acf)- In the limit of F >> Acf the IV state, with N=6, takes over. 

An encompassing phase diagram was proposed (see Fig. [TjD) computing the difference {Tqs — Tcf) versus the sum 
(Tgs +Fcf) of those parameters. According to experimental evidences from the temperature dependence of electrical 
resistivity {p{T, x)), the LM regime correspond to systems with Tgs <Tcf and the ItM one to those with Tgs > Tcf 
[17j . Recently, some non magnetic Ce compounds were reported to follow Coqblin-Schrieffer model predictions for a 
four fold = 4 GS [TM^ . Nevertheless, the limit between a A^ = 4 GS and a quasi-quartet (i.e. two doublets with 
small but comparative Tgs > Tcf) are difficult to be discriminated experimentally. 

II. DIFFERENT TYPES OF MAGNETIC PHASE DIAGRAMS 

With the growing interest on quantum critical phenomena, a number of theoretical models were proposed to 
describe low lying energy excitations related to the Tmo ^ physics, see for example Refs. [S] l2TI - [23] and references 
therein. Since the scope of the models is generally focused on microscopic mechanisms and they are applied on 
specific exemplary cases, there is an absence of an encompassing criteria able to detect or discard new candidates 
and, eventually, to recognize new alternative behaviors. Thermodynamic postulates provide the proper tools for 
such a purpose because of their simplicity and universality. Notably, the constraints imposed by the third law of 
thermodynamics (e.g. the T = singularities in thermodynamic parameters) are some times left aside. 

Another relevant aspect concerns whether there is any condition to be fulfilled at finite temperature for a real 
system to actually reach the Tmo ~^ limit. It is evident, for example, that many Tmo{p) phase boundaries are not 
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FIG. 3. (Color online) Magnetic contribution to specific heat (after respective phonon subtraction) divided temperature of 
three exemplary compounds approaching their respective critical regime, (a) AF-Celns-j^Sui,, (b) FM-CePdi_a:Rhj; and (c) 
AF-CePd2(Gei_a;Sia;)2 after Dash-dot curves indicate the tendency to a nearly constant value of Cm/T maxima within 

the a;o < a; < Xcr region. Notice the logarithmic T dependence. 



properly checked to involve the corresponding entropy of the ordered phase (Smo) despite their Tmo{P) are naively 
extrapolated from T > IK down to e.g. their respective superconductive domes at T^c < 1. Similar questions apply 
to concentration driven systems since the arbitrarily monotonous extrapolations of Tmo{^) exclude any variation of 
dTMo/dx, even in the range where thermal and quantum fluctuations start to compete in energy at low temperature. 

The amount of intermetallic compounds claimed at present to be candidates for quantum critical behavior allows, 
and even requires, such a comparative analysis of their thermodynamic properties in order to distinguish between 
reliable candidates and potential ones. Furthermore, novel alternative behaviors to the quantum phase transitions 
(QPT) might be missed due to the mentioned simplistic extrapolations done in low temperature phase diagrams. 

The simplest starting criterion for such an analysis is to take into account all possible alternative scenarios for 
the decreasing entropy of the magnetically ordered phase (S'a/Oi computed from zero up to Tmo)- There are three 
possible alternatives: 1) Smo —J' as Tmo — >■ 0; 2) Smo at finite temperature due to an eventual transference of 
degrees of freedom to a non-magnetic component, and 3) Smo does not decreases proportionally to Tmo, producing 
a entropy accumulation at low temperature. These three cases are schematically represented in Fig. [2] where each 
case is identified with the respective number. 

The exemplary systems selected for this study are Ce-lattice compounds, where composition or chemical pressure 
effects are produced by alloying Ce-ligand atoms. This criterion preserves the periodicity and local symmetry of the 
Ce net in order to minimize side effects such as disorder [211 ■ Hence, hereafter we will refer to 'alloyed compounds' 
as those where only Ce-ligands are doped or substituted. Otherwise, the few experimental results extracted from a 
Ce diluted system will be explicitly indicated. 
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FIG. 4. (Color online) (a) Example of vanishing specific heat jump at finite temperature as a function of pressure after Ref. |32) 
(open symbols) and Ref. [34] (continuous lines). Dashed curve: non-magnetic HF reference for e.g. 2.4 GPa isobar, (b) Various 
vanishing specific heat jumps as a function of Tm(p)- Notice that Tn{p) decreases to the right a pressure increases. Open 
symbols (left axis) extracted from ac-specific heat and full symbols (right axis) from standard heat pulse techniques, with 
Celus comparing both techniques. 



1. Ce systems with Smo as Tm — >■ 



This is the thermodynamic condition to be fulfilled for reaching a zero temperature QCP. For this group, we 
have selected some exemplary Ce-lattice alloyed compounds whose specific heat were measured down to very low 
temperature, see Fig. [3] The common feature of these experimental results is that the respective specific heat 
jumps {ACm/T) at Tmo first decrease and then broaden levelling off around 2: « xq [23 . This change of regime 
at a: « a;o is associated to a change of slope of the phase boundary at x ~ xq whose implications are discussed in 
the following Section. The Cm/T maximum for x > xq can be analyzed within the Ginzburg-Landau theory for 
second order transitions: Cm/T = where a and h are the coefficients of the free energy expansion G{ip,T) = 

Gq{T) + a{T)4>'^ + b{T)tp'^. The tendency to a constat value of Cmax/T indicates that the G{il', T) dependence on the 
ratio is locked and consequently the entropy evaluated up to that maximum decreases monotonously, following 
a sort of law of corresponding states [26] with the critical point at T = 0. This is a necessary thermodynamical 
condition for any system to reach a QCP. Based on this analysis, one may include into this group the well known 
system CeCue-^Aua: [17] and Ce(Pdi_a;Nia;)2Ge2 [H] because they exhibit the same fatures. one may even propose 
a potential candidate like Ce(Pdi_a;Rha;)In [29] that, to our knowledge, was not yet investigated down to sufficiently 
low temperature. Among the pressure driven systems, one finds CePba [30j showing this type of behavior up to 7GPa 
where a change of magnetic structure occurs. Since this class of magnetic phase diagrams are directly related to 
quantum critical phenomena, we shall discuss in detail some selected experimental results in Section III. 
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2. Systems with Smo — > ai finite temperature 

A second class of magnetic Ce-lattice systems behave quite differently because Smo decreases faster than T^v 
extrapolating to Smo = at finite temperature, being all experimental examples AF. This type of behavior is 
frequently observed in pressure driven stoichiometric compounds and the relevance of this class of phase diagrams 
arises from the fact that many of them exhibit a superconductive GS under pressure |31| . 

The appearance of a superconductive dome is currently related to the extrapolation of the AF phase boundary. 
However, a thermodynamic analysis of those phase boundaries (mostly constructed from resistivity measurements) 
reveal that such a putative extrapolation of Tjv(p) down to Tsc is quite arbitrary. It is known that to extract 
absolute values from heat capacity measurements under pressure exceeding 1.5 GPa is quite difficult. Nevertheless, 
measurements up to 1.2 GPa, e.g. in Celna :32] and CePd2Al2Ga [53], are good reference for those performed at 
higher pressure, e.g. Celua [31], CePd2Si2 [35 and CePd2Ge2 [SB]. These results provide convincing information 
to recognize their distinct behavior respect to those described in the previous subsection (see Fig. |3|. Since very 
high pressure results are given using arbitrary units, a quantitative evaluation of SmoIp) variation is not possible. 
Alternatively one may evaluate the relative decrease of the specific heat jump at Tmo (ACMo(l'))taking into account 
that ACmo — >■ imphes that Smo 0. 

The common feature of these results is that they show a progressive transference of the magnetic degrees of freedom 
to a non-magnetic HF component in the region where Tn{p) decreases. As an example, the Cm{T) results obtained on 
Celna [32l [34] is presented in Fig. |4^. The comparison with the other compounds is done using the relative variation 
of ACmo driven by pressure as depicted in Fig. In that figure (right axis) the quantitative comparison result is 
made using the results from CePd2Al2Ga [55] . 

The main conclusion extracted from Fig. |4jD is that in all these compounds the AC,„(Tjv) jump vanishes at finite 
temperature independently of their eventual superconductive GS, with the phase boundary vanishing some degrees 
of temperature above Tsc- In the case of Celna, this behavior is confirmed by ^^^In-NQR measurements '37' under 
pressure. Detailed electrical resistivity (p) measurements performed on CePd2Si2 |38j also supports this observation 
since the temperature derivative dp/dT (expected to be qualitatively related to the specific heat [38]) shows a jump 
at T = T/v which also vanishes at finite temperature. 

Noteworthy is the fact that these pressure driven superconductors show a similar T/v(0) x po — 2.9 ± 0.2KGPa 
product [39] ■ being Tn{0) the AF transition temperature at ambient pressure and po the pressure where the tran- 
sition vanishes. The origin of this empirical relation is not yet elucidated, but it is confirmed by a number of 
non-superconductive systems which do not verify that product. Among then, the compounds undergoing a maximum 
of their phase boundaries as a function of pressure (e.g. CePd2Al2Ga [33] or CePt [4^) clearly do not present a 
superconductive phase. 

It should be mentioned that the competition between magnetism and superconductivity observed in CeRhlus for 
example [_41j cannot be included within this group because its magnetic transition is of first order and there is a 
coexistence of both phases below T 2K and it merits its own analysis. 

Concerning Ce-lattice alloyed systems, one can mention Ce(Rhi_a;Ru2,)2Si2 [12], Ce(Rhi_2;Pdi^)2Si2 and 
Ce(Rhi_j;Ru2,)3B2 [12] which show equivalent vanishing process of their magnetically ordered degrees of freedom. 
Interestingly, two ways of transfer can be distinguished between ItM (large Fermi surface) and LM (small fermi sur- 
face) magnetic systems. Within the former group (e.g. Ce(Rhi_2,Pd2,)2Si2) the ordered state seems to build up as a 
condensation of degrees of freedom from the 4/ narrow band heavy quasi-particles because an entropy compensation 
is observed respect to a high temperature (T > Tjv) non-Fermi-liquid (NFL) component which does not change with 
concentration |12| . On the contrary, in the second type the transference of degrees of freedom occur between two 
different components because the NFL component increases at the expense of the exhausting magnetically ordered 
ones. Preliminary C,„(T) results on Ce(Coi_2:Fe2,)Si alloys clearly show this transference between two systems be- 
cause AC„i(T/v) vanishes around T w 4.5 K with a coincident rising up of a C„i/T on — Ln(r/To) contribution. No 
superconductive GS is expected in these alloyed systems because of the Ce neighbors random distribution. Neverthe- 
less, the non-magnetic stoichiometric limit {x = 1) of Ce(Rhi_2,Ru2:)3B2 shows low temperature superconductivity 

m- 

3. Systems with critical entropy accumulation as Tjv decreases 

According to thermodynamics, if Smo does not decrease proportionally to Tmo{x) (like in the mentioned first 
class) but with lower ratio, a limT->.o Smo > would occur. In that case, an entropy bottleneck occurs and the 
magnetic phase boundary shall end in a CP at finite temperature where a first order transition drives the system to 
Smo 0. Such a situation is observed in the compounds included in Fig. [5} Ce2(Nii_2,Pdj;)2Sn [H] and URu2Si2 
|45| . The former is a recently studied compound driven by Ce-ligands alloying, whereas the latter is the well known 
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FIG. 5. (Color online) (a) Specific heat divided temperature of the concentration dependent Ce2(Nii_a;Pdi:)2Sn system after 
Ref.[44] and (b) referent field dependent URu2Si2 after Ref.|45] 



0.0 



0.2 



x/x^ [Sn,Pd] 

0.4 0.6 0.8 



1.0 



00 3 

CD 

O 



CO 



Cein, Sn 

3-x X 



,^ Ce2(Pd,Ni)^Sn 




B[T] 



FIG. 6. (Color online) Comparison of the entropy gain Smo up to Tm between two types of behaviors: Celns-ajSn^ of the 1st. 
class and Ce2(Nii_i,Pdj;)2Sn, and URu2Si2 of the 3rd. class. Upper x-axis for concentration x dependent systems, and the 
lower X-axis for the magnetic field dependent one. Notice that Ce2(Nii_2,Pda;)2Sn contains two Ce-at. per formula unit. 
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U compound driven by magnetic field and taken as a referent system for our purpose. It is worth to note that the 
Cm(a;orB)/T variation of the maxima are described by the same function: 6.5 and 7.2/T respectively, both drown 
in Fig. [5] In contrast to the behavior discussed in subsection II- 1 where Cm/Tmax becomes constant, in this case is 
the Cmax{TN) value that remains nearly constant till the first order transition occurs (at _B « 33T in URu2Si2). The 
first order character of the transitions are recognized from the value of the C,„(T;v) maximum which clearly exceed 
the cx l/T function. 

There are also striking coincidences in Smo concerning their similar and nearly constant values that can be ap- 
preciated in Fig. [6] The fact that different control parameters applied on different systems produce the same effects 
can be taken as a fingerprint for the universality of this behavior. For comparison, in Fig. [6] the Smo values obtained 
from CeIn3_a;Sna;, which belongs to the first group, are also included to show that approaching the critical region the 
entropy of these two compounds exceed that of one following the low of corresponding states with a CP at T = 0. 

III. PECULIAR PROPERTIES OF CE-LATTICE SYSTEMS ACCESSING TO Smo ^ 

At the time when the phase diagram presented in Fig. [l] was proposed, no quantum fiuctuation effects were yet 
identified and the usual low temperature limit for magnetic studies sa 1 K didn't provide evidences for such a scenario. 
Thus, the Tmo{x) ^ limit was usually extrapolated following the classical negative curvature to the concentration 
defined as xq in Fig. [T] Later on, lower temperature measurements made evident that approaching xq a change of 
curvature (as presented in Fig. [7| occurs around T^^ sa 2.5 K [47]. Hence, the actual quantum critical concentra- 
tion Xcr does not coincide with the Tmo extrapolation to xq but it occurs at higher values of x. Nevertheless, 
remains a relevant concentration because it characterizes the high temperature region dominated by classical thermal 
fluctuations. To our knowledge, this change of slope was not observed in phase boundaries driven by pressure nor 
by magnetic field. However, there is a striking coincidence in the fact that the phase boundaries driven by pressure 
vanish at similar or higher temperatures than T*^^. 

Apart from the mentioned modification of the Tmo{x) curvature, the change of regime at a;o coincides with other 
unpredicted features occurring around that concentration. Among them, there is a first order transition T/ (x) observed 
the compounds included in Fig. [Sj c.f. AF CeIn3_a;Sna; [IS] and CePd2Ge2_2:Sia; [15] and FM Ce(Pdi_a;Rh3,) single 
crystals [SU]. The exemplary case of CeIn3_a;Sna; is depicted in Fig. [s] and will be discussed in detail in Section 
III-B. Other Ce systems showing a satellite first order transition in similar context are Ce(Pdi_a;Rha;)In [53] and 
Ce(Cui_xNia;)2Ge2 [51]. All these transition shows a similar sharpness, only depending of the quality of the sample. 

A. Magnetic Phase Boundary within the Quantum Critical regime 

The observed change of regime can be explained taking into account the competition between the decreasing energy 
of thermal fluctuations and the temperature independent energy of quantum fluctuations. While the classical transition 
extrapolates to x = xq with coherent thermal fluctuations of the order parameter decreasing with temperature, the 
critical fluctuation associated to a QPT respond to a quantum-statistical description [3^. As mentioned before, only 
composition driven phase diagrams show this clear change of regime at T'^^K (see Fig. [t^) . In composition driven 
systems, the formation of static disorder fluctuations (or "rare regions" [TI) was proposed to explain the shift of the 
critical regime from xq to Xcr- 

The occurrence of "rare regions" together with the so-called Griffiths effects [3] are mostly addressed to describe 
composition driven systems because of their intrinsic possibility of local disorder attributed to all alloyed systems. 
As it was discussed in Subsection I-C, there is a clear difference between structural and chemical potential pressure 
effects because the former implies a random (i.e. disordered) Ce-ligand atomic volume distribution. This effect 
cannot be nai'vly extrapolated to chemical potential variation introduced by neighbor elements with nearly equal 
atomic volume substitution. This effect is therefore better described as an electronic topological distribution. In fact, 
the change of regime occurs: i) quite suddenly and only close to xq, ii) in different compounds with quite different 
relative concentrations and iii) without a further broadening at higher concentrations (52) . Moreover, the first order 
transitions at Tj (see e.g. inset in Fig. [s]), show similar sharpness in all cases despite of their different concentration 
regions excluding any atomic disorder as a relevant factor. 

All these experimental evidences confirm that there is a crossover between two distinct regimes at similar tem- 
perature. Beyond that region, QPT mechanisms dominate the low energy scenario producing drastic changes in the 
nature of the magnetic phase boundary as it will be analyzed in detail in the following subsection. Some of those 
effects can be observed in Fig. [Tja for two AF compounds where Tm o<:\ x — Xcr \ and in Fig. [?}; for a FM one, where 
Tc{x) decreases asymptotically till it collapses to zero at very low temperature. Notice that for the x > Xq region the 
phase boundary is label as Tm instead of Tmo in order to distinguish them as belonging to different regimes. If one 
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FIG. 7. (Color online) (a) Magnetic phase diagrams of three exemplary compounds covering an extended range of temperature, 
(b) and (c) show the detail of the quantum critical region for: (b) two AF compounds with linear Tm{x) dependence and (c) 
with an asymptotic phase boundary of a FM compound [IS] (notice the logarithmic Tc axis). 



describes these magnetic phase diagrams within the pattern proposed in Fig. [T] one may recognize that the systems 
included in Fig. [7)d belong to the xq < xk class whereas that from Fig. ^ to the xq > xk ones. 

To gain insight into this unexplored range of concentration between Xq < x < Xcr, we will address our attention to 
the phase diagram and the related thermal properties of Celna-ajSn^;. In Fig.|8]we shown a detail of the phase diagram 
around the critical region where three phase boundaries converge: i) the classical AF-Tjv(a;) dominated by thermal 
fluctuations, ii) Tm{x) dominated by quantum excitations, and iii) the dome of a first order transition Ti(x) curve. 
The different nature of Tm respect to the classical AF-transition at Tpf(x < xq can be clearly appreciated in the inset 
of Fig. [sjfrom the flattening of the transition by magnetic field (up to i? = 4T) measured on sample x = 0.45. The 
same behavior is observed at higher concentration down to the milikelvin temperature range [54 . As a comparison, 
one can mention that magnetic field applied to an x = 0.25 < sample [55] shows that a very high field (about 40 T) 
has to be applied to flatten the Tjv transition, which then decreases in temperature before to vanish. Coming back 
to the studied x = 0.45 sample, the jump at Tm can only be observed analyzing the derivative of Cm{T) as shown 
in the lower part of that inset. To our knowledge, the CR scenario presented in Fig. |8] was only recently pointed 
out by theory for itinerant FM systems [56|. However, our observations are also include quite localized AF systems. 
Although a first order dome occurs close to CR similarly as proposed by theory, the following Tm{x) transitions 
resamble Pippard's third order ones [57) . 
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FIG. 8. (Color online) Detailed phase diagram of CeIn3_a;Sna; around the critical region CR. One can identify different phase 
boundaries: second order transitions Tjv for x < 0.3, a first order dome Tj around CR (0.3 < a: < 0.45) and the Tm boundary 
within the pre-critical region {xo < x < Xcr)- Inset: Magnetic field effect on Tm and Ti transitions on sample CeIn2.55Sno.45 
after Ref.[3S]. Right axis (lower curves) : suppression of the Cm{TM) derivative jump for _B = 4T. 



B. Thermodynamic behavior of Celna-ajSria; within the xq < x < Xcr range 

In this subsection we will analyze the Cm{T) and entropy {Sm{T)) contributions of the Celna^ajSn^: system, which 
was investigated in detail within the xq < x < Xcr range |48j . Taking profit of the linear variation of Tj\/ with one 
can compare the C'm{T)/T dependence of different samples by normalizing the temperature as AT(x) — T~Tm{x). In 
Fig.[9]we show that comparison performed on seven samples covering the concentration range between 0.41 < a; < 0.80. 
Notice that with this definition of AT, the magnetically ordered phase corresponds to the negative range of that 
parameter, see the upper x-axis of Fig. |9] There, one can see how the Cm{T > Tm)/T tails of the alloys belonging to 
the 0.45 < X < 0.70 range converge into a unique curve. In order to remark the validity of this scaling, we have also 
included in the figure the results obtained from x = 0.41 and 0.80 samples, placed beyond the limits of the pre-critical 
region, which clearly deviate from the scaled ones. 

In Fig. [9] the vertical line at AT = line splits the Cm/T contribution into two parts, the T < Tm one (hereafter 
label as Cm/T) and the tail at AT > hereafter identified as Cnfl/T because of its NFL behavior (c.f. oc Ln{T /Tq)). 
Notably, also the Cm/T contribution for the samples within this pre-critical region overlap each other at AT < 0. 
Samples x = 0.41 and 0.45 show the weak peak related to a first order transition discussed in the previous subsection. 
The relevant feature is that the Cm/T overlap allows an extrapolation of Cm/T — ^ to AT w — 4K which is also 
independent of concentration. We remind that a AT < value does not correspond to a negative temperature but 
simply to a common extrapolation to a zero value of the C^/T contribution. 

The key parameter to describe this peculiar behavior of the specific heat is its associated entropy, evaluated as 
Sm — J Cm/TdT. According to the definition proposed for AT, one may split the total measured values as Sm = 
Smo + Snfl, being Sm the contribution of the T < Tm phase (i.e. AT < 0) and Snfl the one from the NFL 
tail for AT > 0. For the following analysis we take as reference the entropy variation of sample x — 0.41 because 
it contains largest Sm contribution among the samples included in Fig. [9] As it can be appreciated in the figure, 
the Sm{x = 0.41) contribution slightly exceeds 0.2RLn2 whereas Snfl reaches ~ 0.6RLn2 (see inner right axis). 
Noteworthy, the full RLn2 value is only reached once the extrapolation to the Cm/T = value at AT w — 4K is 
included, as depicted using the lowest 'AT -I- 4' and outer ^Sm + Snfl right axes in Fig. |9] Since Snfl ~ 0.6RLn2 
does not change with concentration, but Sm — as a; — >■ Xcr, one concludes that about 40% of the RLn2 entropy is 



missed as Tm ^ 0. that difference is illustrated in Fig. 10a, by comparing the S„i{T) variation from samples x — 0.15 
(with full entropy) and x = 0.70 ~cr- 

It is evident from Fig. |9]that the decrease of SMix — > Xcr) does not imply a transference of entropy to the NFL 
phase because Snfl is independent of concentration within this concentration range. Consequently, those degrees 
of freedom are missed without any change of SnfLi which remains unchanged with a 60% of the RLn2 value. We 
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FIG. 9. (Color online) Overlap of Cm/T'(0.45 < x < 0.70) curves plotted versus a shifted temperature AT = T — Tm (upper 
and left axes) after Ref.|25|. The dashed curve represents the concentration independent entropy gain for 0.45 < a; < 0.70 
samples, discriminated between the T < Tm range Sm and the paramagnetic one Snfl (low T and inner S axes). Full RLn2 
entropy is computed from AT — — 4K (lowest 'AT + 4' and outer 'Sm + Snfl' right axes), see the text. 

resume this situation in Fig. [T0| 3, where measured Sm up to T = 20K and 5*0 = RLn2 — Sm are represented as a 
function of concentration. At the QCP, Sm — > Snfl because Sm — >■ and Sq — >■ iO%RLn2. 

Simplistic explanations looking for a some extra entropy contribution at higher temperature fail because it would 
imply a discontinuous transference of entropy from the T < Tm to temperatures above 20 K according to Fig. |9] We 
recall that in Celna the crystal-field excited quartet lies at high enough energy (« 100 K [58]) which guaranties no 
contribution to this analysis. Neither a Kondo temperature increase can be argued because the temperature of the 
maximum of the electrical resistivity (T^oa;) remains unchanged between xq and Xcr at T^^a; ~ 19K |48j . 

The lack of entropy showed by this Ce compound is not an exception because in the cases where this type of 
analysis was performed the RLn2 value for a doublet GS was never reached. For example, the compounds showing 
a Cm/T (X Ln{T/To) dependence only reach a 0.54i?Ln2 value ^59j. Since the mentioned Ln{T/To) dependence 
corresponds to one of the possible scenarios for QCPs predicted by theory [3] , one infers that this deficit in the entropy 
or the consequent arising of remnant entropy at T = (Sq) is intrinsic to the NFL phenomenology approaching that 
point. 

The question arises whether the low temperature x — >■ Xcr behavior is governed by low dimensional fluctuations as 
it was observed in CeCus+j^Aui-j; |60j . Unfortunately neutron scattering studies on Celna-ajSn^; are not available 
because of the strong neutron absorption of In nuclei. Alternatively, one may check whether any evidence for low 
dimensional fluctuations can be extracted from thermodynamical results. For that purpose one can evaluate the 
internal energy, Um{x) and Sm{x) for Tm < T < oo in samples around that concentration, and compare them with 
Ising and Heisenberg model predictions 61j for 1, 2 and 3 dimensional systems with different lattice structures (i.e. 
coordination number). Such analysis showed that the those thermodynamic parameters nicely fit into the predicted 
values for a 2D- Ising quadratic layer [52] , 



C. Thermal Expansion 

In order to confirm that the anomalous evolution of the entropy approaching the critical point as due to an 
intrinsic effect, a complementary thermodynamic parameter sensitive to this phenomenon has to be investigated. 
Such alternative is provided by thermal expansion f5(T,x) measurements because they are related to the entropy 
through the Maxwell relation —dS/dP — dV/dT — /3.Thus an anomalous Si^{x — x^r) dependence should have a 
replica in Vo(a; — Xcr) as T ^ 0. In this case, the effective pressure is originated in the chemical pressure produced 
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FIG. 10. (Color online) a) Comparison of the temperature dependence of Sm between the alloy x — 0.15 belonging to 
the classical region and one lying on top of the critical concentration x = 0.70 showing a deficit of ~ 40%. b) Left axis: 
concentration dependence of the entropy Sm = Sm + Snfl of Celns-^^Sni, samples measured up to 20K, and the zero point 
entropy So computed as the difference respect the total expected value RLn2. Xcr highlights the critical concentration. Right 
axis; relative change of volume at T — > for x > 0.4 samples, see text. 



by alloying. 

The /3(T — )■ 0, x) variation of CeIn3_2;Snj; was studied down to the mK range in the vicinity of the critical concen- 
tration |63j . In Fig. 10 3 we have included the volume variation as Vo{x) for T — extracted from the temperature 
variation V(T) = jpdT. The obtained values were normalized well above any quantum fluctuation effect, i.e. 
AK < T < 8 K taking as reference the x ~ 0.8 alloy which lies beyond the critical point. Both abnormal So{x) and 
Vo{x) dependencies are compared in Fig. 10 3 around Xcr- 



IV. THERMODYNAMIC IMPLICATIONS OF THE T ^ PHYSICS 
A. Third law of Thermodynamics and Remnant Entropy 

Thermodynamic postulates |26j state that the entropy decreases as T — > reaching a finite value 5*0, which is 
not necessarily zero [57] . The implication of this postulate can be viewed in a different way by considering the 
thermodynamic definition of temperature [55], i.e. that T = is reached once dU/dS — 0. Thus, the thermodynamic 
condition for zero temperature corresponds to a zero variation of the internal energy U and not to the value of the 
entropy itself. The Sq = value corresponds to a singlet GS without any other accessible degree of freedom [M] . 
Meta-stable states may eventually decay into such a GS in infinite time (like e.g. amorphous or other structurally 
disordered systems), however frustrated systems or those dominated by quantum fluctuations (like those involved in 
the present study) escape to this possibility. 

The way to prove that So ^ was applied in Section III-C by taking as reference the Sm = RLn2 value associated 
to the magnetic doublet-GS of Celn^-xSnx at T sa 20 K. Interestingly, the absolute reference for the the entropy value 
is taken form high enough temperature where both levels of the doublet GS are equally occupied (i.e. Sm — RLn2). 

A complementary aspect regarding the application of the third law of thermodynamics implies the lim7'_>.o dS/dT. 
The dS/dT = possibility at T = corresponds to the previously mentioned case of a singlet GS, realized in a long 
range ordered state and associated to a positive 5,„(r) curvature (i.e. d'^Sm/dT^ > 0). The simplest examples for a 
dS/dT 7^ (or d'^Sm/dT'^ = 0) is provided by metallic systems whose conduction electrons are described by standard 
Fermi-liquid behavior with dS/dT = 7, c.f. the Sommerfeld coefficient. Heavy fermion (HF) systems simply increases 
the dS/dT slope without changing their physical properties. 

Since a negative Sm{T — )• 0) curvature (i.e. d'^Sm/dT'^ < 0) is not possible because it would imply a singularity 
at T = 0, the question arises whether there is an upper limit for the dS/dT = Cm/T slope or, equivalently, for the 
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FIG. 11. Examples of maximum slope of Sm{T) at T — >■ in real systems, (a) CePdaB [68], CeCue-a^Aux [271 ? ], 
(Cei_yLay)TiGe [SS] and CeAge-zAuz [BS]. (b) Detail at the mK range including: GeCus.gAuo.i (dashed line), CeGu5.8Auo.2 
at different pressures and, for comparison, the fit for CePdaB from T > 0.5 K (dot dashed line). 

7 value in non ordered HF. Such a question is related to the divergence of thermal parameters like specific heat or 
thermal expansion because even a logarithmic divergence at T = implies an infinite slope of dS/dT{T — ^ 0). This 
scenario was proposed by theoretical models [3] and claimed to correspond to experimental results [2 dismissing some 
thermodynamic postulates. 

We will analyze now the existence of an eventual upper limit for dS/dT(T — > 0) and, in the following subsection, the 
thermodynamic consequences of the proposed divergencies at T — >■ 0. In Fig.[TT]we have collected the low temperature 
Sm{T) dependencies extracted from a number of Ce system showing the highest values of Cm/T > 3 J/molK^ for a 
doublet GS {N — 2) independently that they order or not. This comparison includes the well known CeCu6-a;Aua; 
[27] and CeCug-xAga; jM] Ce-lattices. Also the diluted (Ceo.iLao.9)TiGe [66] is included to confirm that this limit 
is related to thermodynamic properties independently of a lattice configuration. Other Ce diluted systems, like 
(Cei_2;La2:)Pt3Si 67J, also exhibits similar high values for 0.8 < a; < 0.95. 

As it can be seen, all experimental results show a sort of envelope curve described by the Sm{T) depen- 
dence of CePdaB [68]. To our knowledge, the record of low temperature entropy was extracted from the Ce- 
diluted (Ceo.o3Lao.97)B6 jSS] which exceeds the low temperature Sm{T) values of CePdaB even after normalized 
by RLnA (remind that CeBg has a N=4 GS). Among the Ce-lattice systems, CeNi9Ge4 shows a very high value: 
Cm/T{T — >• 0) = 5.2 J/molK^ |2D], because its GS is composed by the contribution of two hybridized doublets with 
a Kondo temperature equivalent to the extremely small CP splitting (« 10 K). Also in this case the corresponding 
normalization entropy is RLnA. 

From these experimental evidences we conclude that there is an upper limit for dS/dT{T — > 0) w 4J/molK^ for 
N — 2 GS in Ce systems and consequently a lower limit for the Kondo temperature evaluation based on the thermal 
dependence of the entropy. The question arises whether this is a physical limit for Tk or simply a limit for the 
application of theoretical models. 



Two types of divergencies for thermodynamic parameters were proposed to account for the low energy excitations 
in NFL systems [5], one described by a power law like Cm/T (x and the other by a logarithmic dependence 
like Cm/T oc —ln{T/To). Details concerning the range of applicability of different models exceed the scope of this 
phenomenological work and can be found in many interesting review articles [U [2TH23I \60\ [70] . 

Most of Ce systems obeying a Cm[T)/T power law dependence are described by non fractional values of the 
exponent g, with scarce coincidence with theoretical predictions devoted to systems with short range order parameter 
fluctuations |23[ 170] . Divergent power laws imply a non analytical singularity at T = which is in conflict with 



B. Divergencies at T — )> 
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the third law of thermodynamics. In agreement with thermodynamics, actual specific heat results show a systematic 
tendency to saturation of Cm{T)/T = at the low temperature limit, which is well described by an heuristic 
modified power law jj- — G/{T'^ + A) [5S]. In that formula A represents an energy scale below which the third law 
constraint of Cm with T [25^ becomes dominant. That parameter also allows to evaluate the temperature 
of the crossover between NFL and FL regimes since the later does not diverge at T — 0. The T = limit of that 
formula is obtained computing the limit of 7t^o = G/A. In the case of the CePdsB compound, proposed in Fig. 11 as 
the phenomenological envelope curve, the fitting function obtained between 0.5 and 4K, is C^/T = A.5/{T'^-^ + 1.1) 
which results in a 7t-*.o ~ 4 J/molK^ value. 

The other usual divergence observed in NFL systems has a logarithmic character. Besides the fact that this type of 
algebraic divergence does not imply a singularity at T = 0, even the system showing the highest Cm/T values observed 
at low temperature in CeCus.gAuo.i [50] (included in Fig. 11 1 does not exceed the S'm(T) curve represented by CePdsB 
significantly. The comparative study performed normalizing the temperature with the respective energy scale Tq as 
t = T/Tq ^ shows that a universal function Cm/t = —DLog{t) + Eq x Tq describes all analyzed compounds, with 
D = 7.2 J/molK and Eq accounting for any eventual high temperature contribution (which is zero for CeCus.gAuo.i). 
Once subtracted the non logarithmic contribution Eo x Tq the computed entropy does not exceed 60% of the RLn2 
value. 

Interestingly, there is a common feature in systems showing a Cm/T cx —Ln{T/To) divergency before to reach 
the critical point, that is the specific heat jump transformed into a kink. Such an anomaly is well illustrated by 
CeCu5.8Auo.2 and CeIn2.55Sno.45 samples [53], which requires of the temperature derivative of Cm to observe a a 
discontinuity at Tm as shown in the inset of Fig. |8] That feature is also observed in magnetic field driven phase 
boundary of Sr3Ru207 [71^. Based on the fact that in CeCug-a^Au^: the low energy excitations were recognized from 
neutron scattering measurement to arise from magnetic fluctuations with an effective dimensionality smaller than three 
[72j , one may infer that the mentioned cusp reflects the low dimensionality of the order parameter in this quantum 
critical region. Analyzing the Cm{T) dependence around that transition, its symmetry respect to Tm reminds those 
observed in other Ce compounds in similar conditions, i.e. with Tm < 2.5K j47 . 

Thermal expansion also provides phenomenological information for this scenario from measurements performed on 
Ce(Pdi_i:Cu^)2Si2 in two perpendicular crystalline directions: Uc and aab |73j . In this case one observes that, once the 
ACm(rAr) jump (characteristic for a 3-D mean fleld like transition) is smeared by the magnetic decoupling between 
Ce planes, the corresponding a{T) discontinuity transforms into a cusp before to vanish as a function of x. The 
outstanding aspect of this change is that it occurs at different concentrations depending on the measured direction. 
While in 'c' direction the evolution goes hand by hand with the speciflc heat jump vanishing &t x — 0.2, along the 
'ab' plane it holds up to a; = 0.3 but transforming into a cusp at T « 1 K as expected for a lower dimensionality order 
parameter, see Fig. 10 in Ref . [T^ . 



C. Consequences on the lower Tk limit determination 



Further consequences arise from thermodynamical constraints imposed by the third law concerning the evaluation 
of the Kondo temperature from the entropy, e.g. 5„i(Tr-) — 2/3i?Ln2 [TH. If there is an upper limit for the dSm/dT 
derivative, a consequent lower limit occurs on the value oiTu extracted from SmiT). From the envelope curve proposed 
in Fig. |11| such a limit would be w 1.3 K. Similar situation occurs with models extracting Tr- from the speciflc heat 
jump at Tjv as Tk/Tm oc ACo/^Cm [ZS], being ACq = 1.5i? the reference value from mean field calculation of ACm 
for a doublet GS. In this procedure the fixed 1.5R value contradicts the constraint imposed by the mentioned law of 
corresponding states which requires that ACm — >■ as Tjv — >• affecting the application of the ACq /ACm ratio for 
a Tr- evaluation. Another current criterion to evaluate the Kondo temperature is to compute Tx cc I/70 [15 . Also 
in this procedure the minimum value of Tk is limited by the empirical maximum of liuiT^o dSm/dT — 70 observed 
in Ce systems. Whether this low Tx limit arising from thermodynamic conditions on 70 and ACm is intrinsic to 
the Kondo effect itself is an open question. Theoretical models remark the possibility of a quenching of Kondo effect 
approaching a QCP [1TJ[21], but not based on thermodynamic constraints. In any case, these considerations warn on 
the accurate application of theoretical models not accounting for thermodynamical constraints on real systems. 



V. CONCLUSIONS 



This comparative analysis of the low temperature properties of Ce-magnetic systems show the power of thermody- 
namic parameters in recognizing different types behaviors, in particular the not fully profited information extracted 
from the entropy at T — > 0. Moreover, the third law of thermodynamic provides universal sine quibus non conditions 
for real systems to approach a zero temperature QCP, independently of a priori model hypothesis. It is important to 
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distinguish between candidates to present QCP and other with ahcrnativc behaviors because new relevant physical 
phenomena might be missed due to a misleading low temperature extrapolation. 

A significant amount of experimental evidences were analyzed at the light of these criteria, which allow to conclude 
that at least three types of phase diagrams can be clearly distinguished. Depending on the behavior of the Tmo phase 
boundaries, those phase diagrams can be sorted as follows: i) those where the phase transition fulfills the conditions 
to be driven to T = 0, ii) those whose phase boundaries vanish at finite temperature because their MO degrees of 
freedom are progressively transferred to a non magnetic component, and iii) those ending in a critical point at finite 
temperature. 

In the first case the possibility to reach a QCP is supported by the continuous decrease of the Smo entropy, which 
extrapolates to zero as Tmo 0. Despite of its monotonous decrease, the phase boundary driven by alloying Ce- 
ligands shows a change of curvature at x = a;o- This behavior is attributed to a change from classical to quantum 
critical of regime since beyond that concentration quantum fliictuations seem to dominate the scenario. Strikingly, 
such a change occurs at similar temperature energy T*-^-^ « 2.2 ± 0.3 K in all studied systems, and below that 
temperature a tendency to saturation of the maximum of Cm{T)/T arises as a distinctive characteristic according to 
a law of corresponding states for a T = critical point. 

A number of distinctive properties were highlighted by a detailed analysis of the thermal properties of the exemplary 
system Cclna-jjSnj;: i) there is an anomalous reduction of about 40% of the entropy respect to reference value i?ln2 
expected for a doublet GS. This missed entropy can be regarded as a zero temperature remanent entropy, ii) it can 
be quantitatively demonstrated that, contrary to current suppositions, the reduction of Smo as Tmo — >■ is not 
transferred to the paramagnetic phase, and iii) at the critical region defined by x = xq and T = T*^^ there is a 
systematic presence of a first order transition and beyond that point the phase boundary Tm{x) changes its nature, 
manifested in a strong dependence on magnetic field and a Cm{T) jump only observed in its temperature derivative. 

No evidences for T — > divergencies arc observed in real systems, instead a progressive saturation of Cm/T is 
observed in those cases described by a power law a finite temperature. Neither those systems with oc —Ln{T/o) 
dependence exceed the empirical upper limit of 70 « 4 J/molK^ for Ce systems with doublet GS. 

The second type of Smo{Tmo) behavior is currently observed in pressure driven phase boimdaries. In this case, 
specific heat results indicate that the phase boundary itself vanishes because of a progressive transference of degrees 
of freedom to the non-magnetic; c;omponent occurring at T > 2 K. Despite of the formation of a superconductive phase 
their magnetic phase boundaries do not reach that transition because it occurs below the 2K threshold. This type 
of behavior is also observed in Ce-ligand alloyed system, but there the occurrence of superconductivity is unlikely 
because of allying effects. 

The distinct characteristic of the third class of phase diagrams is given by the fact that the Cto(Tjv) maxima values 
are found to be constant instead of Cm/T like in the first group. In this case the entropy accumulation as Tjv decreases 
makes the phase boundary to end at a finite temperature critical point. There, a first order transition drops Smo 
down to 0. This peculiar scenario was detected in a system driven by Ce-ligand composition and confirmed by a well 
know U compound driven by magnetic field. Notably both systems coincide in their Smo values. 

To our knowledge, many of these experimental observations were not predicted by current models focused into the 
physics of QCPs. This is probably due to the difficulty of a generic treatment of thermodynamic parameters like 
entropy or the specific heat jump in a region dominated by a complex spectra of quantum excitations. 
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